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We study the mechanism of the mass enhancement in an intermediate- valent regime of heavy-fermion materials. We 
find that the crossovers between the Kondo, intermediate valent, and almost empty /-electron regimes become sharp 
with the Coulomb interaction between the conduction and / electrons. In the intermediate-valent regime, we find a 
substantial mass enhancement, which is not expected in previous theories. Our theory may be relevant to the observed 
nonmonotonic variation in the effective mass under pressure in CeCu2Si2 and the mass enhancement in the intermediate- 
valent compounds a-YbAlB4 and y8-YbAlB4. 
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The heavy-fermion phenomenon is one of the most remark- 
able consequences of a strong electron correlation. In some 
heavy-fermion materials, the effective mass of electrons be- 
comes a thousand times as large as the free-electron mass. 
Such heavy electron mass is due to the renormalization effect 
on the hybridization band by the strong Coulomb interaction 
U between localized /-electrons. 

After the discovery of the superconductivity in the heavy- 
fermion compound CeCu2Si2,'' several heavy-fermion super- 
conductors have been investigated. Since the onsite Coulomb 
interaction is strong in a heavy-fermion system, superconduc- 
tivity is expected to be unconventional, i.e., other than the s- 
wave, and has been one of the central issues in the research 
field of solid state physics. In many cases, superconductivity 
takes place around a magnetic quantum critical point, where 
the magnetic transition temperature becomes absolute zero. 
Thus, the superconducting pairing interaction is supposed to 
be mediated by magnetic fluctuations in these systems. 

However, in CeCu2Si2^^ and CeCu2Ge2,^^' superconducting 
transition temperatures become maximum in high-pressure 
regions far away from the magnetic quantum critical points. 
In addition, the superconducting region splits into two regions 
in CeCu2Sii.8Ge() 2.'*^ Thus, the superconductivity in the high- 
pressure region in these compounds is difficult to be under- 
stood by the magnetic fluctuation scenario, and the supercon- 
ductivity mediated by valence fluctuations is proposed.^'''* In 
these compounds, the effective mass, deduced from specific 
heat measurements or the temperature dependence of electri- 
cal resistivity, decreases rapidly at approximately the pressure 
where the superconducting transition temperature becomes 
maximum.^' The effective mass m* in heavy-electron sys- 
tems is closely related to the valence of / ions:^' 

m* l-tif/l 

— = -, — (1) 

m 1 ~ «/ 

where m is the free-electron mass and n f is the number of / 
electrons per site. This relation is derived for the periodic An- 
derson model (PAM) with t/ ^ oo by the Gutzwiller method. 
Thus, m* decreases as n/ decreases. In Ce compounds, n/ de- 
creases under pressure, since the /-electron level e/ in a pos- 
itively charged Ce ion surrounded by negatively charged ions 



becomes higher and also the hybridization matrix element V 
increases. Therefore, we expect that a sharp change in or 
large valence fluctuations play important roles in the super- 
conductivity in these materials. 

However, eq. (1) is derived for the ordinary PAM, which 
does not show a sharp valence change. Moreover, the effec- 
tive mass has a peak in CeCu2Si2 under pressure before the 
superconducting transition temperature becomes maximum.^' 
Such a nonmonotonic variation in the effective mass cannot 
be expected from eq. (1). Note also that, in CeCu2Ge2, the ef- 
fective mass shows a shoulder structure before superconduct- 
ing transition temperature becomes maximum.'^ This shoul- 
der structure may also become a peak if we can subtract the 
contributions of magnetic fluctuations, which are large in the 
low-pressure region. These peak structures may be explained 
by a combined effect of valence fluctuations and the renormal- 
ization described by eq. (1),^-' but the applicability of eq. (1) to 
a model with large valence fluctuations is not justified. Thus, 
we should extend eq. (1) to a model that shows a sharp valence 
change to understand the superconductivity in CeCu2Si2 and 
CeCu2Ge2 coherently by the valence fluctuation scenario. 

Another important recent issue on the heavy-fermion phe- 
nomenon is the heavy-fermion behavior in the intermediate- 
valent compounds Q'-YbAlB4 and yS-YbAlB4."> yS-YbAlB4 is 
reported to show superconductivity at a very low tempera- 
ture.'^' Although both compounds show heavy-fermion be- 
havior, the valences of Yb ions are +2.73 for Q'-YbAlB4 and 
+2.75 for y8-YbAlB4.'^> Thus, the hole numbers in the / level 
are «/ = 0.73 and 0.75 for ci;-YbAlB4 and yS-YbAlB4, respec- 
tively. With such Hf 1, heavy-fermion behavior is not ex- 
pected from eq. (1). 

In this research, we study an extended periodic Anderson 
model (EPAM) with the Coulomb interaction Ucf between the 
conduction and / electrons, which induces sharp valence tran- 
sitions, by the Gutzwiller method. We extend the Gutzwiller 
method for the PAM developed by Fazekas and Brandow'"' 
to the present model. This extension is straightforward but the 
formulation is lengthy, and here we show only the obtained re- 
sults. The details of the derivation will be reported elsewhere. 
Although the EPAM has been investigated by some numeri- 
cal methods in recent years, ''*"'^' the effect of Ucf on the mass 



1 



J. Phys. Soc. Jpn. 



LETTERS 



enhancement is not yet clarified well. 
The EPAM is given by'^' 

ko- icr kcr 

(2) 

i i(T(r' 

where Cka- and fka- are the annihilation operators of the con- 
duction and / electrons, respectively, with the momentum k 
and the spin cr. «cio- and n/,o- are the number operators at site 
i with cr of the conduction and / electrons, respectively, ek is 
the kinetic energy of the conduction electron. In the following, 
we set the energy level of the conduction band as the origin 
of energy, i.e., 2* = 0. We set U —> oo, since the Coulomb 
interaction between well-locaUzed / electrons is large. 

We consider the variational wave function given by It/') = 
PffPcfl^), where P/f = Ylil^ ~ n/i^nfn] excludes the dou- 
ble occupancy of / electrons at the same site, and Pc/ = 
Ylia-cr'i^ - (1 - g)nci,Tnfia-'] is iutroduccd to deal with the on- 
site correlation between conduction and / electrons.'"*^ g is a 
variational parameter. The one-electron part of the wave func- 
tion is given by |0> = Uk<kf,Acl^ + aik)flj\0}, where kf is 
the Fermi momentum, |0) denotes vacuum, and a{k) is deter- 
mined variationally. Here, we have assumed that the number 
of electrons n per site is smaller than 2. 

Then, we apply Gutzwiller approximation. Here, we intro- 
duce the quantity d^o- = YjiiridAn fn + nfii))/L, where (■ ■ ■ ) 
denotes the expectation value and L is the number of lat- 
tice sites. In evaluating expectation values by Gutzwiller ap- 
proximation, we determine dar, which has the largest weight 
in summations. The result is = [<ica-(l ~ «/ ~ "eu + 
dca-)]/l(nf - dca-){nccT - dca-)], where rica- = Y^Mda^lL and 
Uf - Yja^fa- - llia-('^fia-)/L. This is the same form as that in 
the Hubbard model, '^^ if we regard rica- as n^, «/ as n^, and 
dca- as d^, where and d^ are the numbers of cr-spin elec- 
trons and doubly occupied sites per lattice site, respectively, 
in the Hubbard model, and & denotes the opposite spin of cr. 

In the following, we assume a paramagnetic state, i.e., 
rifa- - nf/2, rica- = tic 12 - in - n/)/2, and dca- - djl, and 
optimize the wave function so that it has the lowest energy. 
In the following, we regard c/ as a variational parameter in- 
stead of g as is done in ordinary Gutzwiller approximation. 
Then we find that a{k) = 2Vi/{ef -h + [(e/ - h)^ + W^] 
where V2 = ^|qV\ = X qcjV and Ik - qc^k- ?/ is the 
renormalized /-level obtained by solving integral equations, 
as we will show later. The renormalization factors are given 
hyq = \n]{nc-d){l-ncl2){l-nf-ncl2+dl2)]l[{l-nfl2){\- 
nf)nc{nf-d/2)^],qcf - {l+{d(nc-d)/[{nf-d/2){l-nf-nc/2+ 
d/2)]}''^/2)inf-d/2)^/[njil-nc/2)laRdqc = qc^ = 

dca-)(l-nf-nca+dccr)V'^ + [dcJnf-dccr)V'^f/[ncJl-nca-)]. 
qctr has the same form as the renormalization factor q^ in the 
Hubbard model'^^ as for the Gutzwiller parameter g. 

To determine nj. If, and d, we solve the following inte- 
gral equations. Uf = n/2 + 1^, €/-€/ = -2V2hqdq^^ /drif - 
(h - h - hef)q~^dqcldnf + AVjhq'jdqcf/dnf, and 
Ucf = -2Vlhqdq-ddd - di - h - h~€f)q-^dqcldd + 
W^hq^j-dqcf/dd. The integrals are given by /i = T^kKkf h/L, 
and // - Zk<kMh - ~eft^im - Iff + AVl\"^}IL for / = 2- 
4. The total energy per site is h + UfCf + (n/2 - - h- 



AVlh + Ucfd. 

We can evaluate expectation values of physical quantities in 
the optimized wave function. Here, we consider the jump in 
the electron distribution at the Fermi level; the inverse of the 
jump corresponds to the mass enhancement factor. The jump 
in ndk) = {c^j^^Cka-) at the Fermi level is given by Ancfe) = 
qcl[\ + qaHkp)]- The jump in nf(k) = {fl^fko) is given by 
■^/(^f) = ^/^fl^(A:F)/[l+<7fl^(A:F)]. The renormalization factor 
qf for an / electron is given by qf - q^p\^^^\l - nf)/(l - 
nf/2), where = {[(«/ - dca-){l - Uf - Hco- + dca-)]^'^ + 
\dca{nc(T - dca-)]^^^}^ l\nf{\ - «/)]. (^^'^^ has the same form 
as q^ in the Hubbard model, '^^ if we regard n/ as n^, Uca 
as n^, and dc^ as d^. In the following, we call 1/An(/:F) = 
l/[Anc(fcF) + An/Cfcp)] the mass enhancement factor. 

Before presenting our calculated results, here we consider 
three extreme cases in the model. First, we consider a case 
with a positively large e/. In this case, «/ ^ and the energy 
is almost the same as the kinetic energy of the free conduc- 
tion band with «<; = n. Second, we consider a case with a 
negative e/ with a large magnitude. In this case, «/ - 1 and 
He - n- \. The energy is approximately given by the sum of 
L[e/+(n-l) Ucf \ and the kinetic energy of the free conduction 
band with Hc = n-\. We call this regime the Kondo regime. 
From the form of the renormalization factors and a(k), the 
mass enhancement factor becomes large as — > 1, which 
is consistent with the previous result on the PAM. Third, we 
consider a case with an intermediate ej with a large Ucf - In 
this case, the / and conduction electrons tend to avoid each 
other, andthusn/-i-ne/2 ^ 1 anAd ^ 0. That is, «/ ^ 2-nand 
Uc -2n-2. Here, we call this regime the intermediate-valent 
regime. In this case, both the / and conduction electrons are 
almost localized, and the energy is approximately L{2 - ri)ef. 
In this intermediate-valent regime, the mass enhancement fac- 
tor becomes large as «/ -h ncl2 1 and d ^ 0. This mass 
enhancement in the intermediate-valent regime is not reaUzed 
in the ordinary PAM and is a result of the eff'ect of Ucf- 

In the following, we consider a simple model of the ki- 
netic energy: the density of states per spin is given by p(e) = 
1 /(2W) fov-W <£<W; otherwise, p(e) = 0. 

Now, we show our calculated results. Figure 1(a) shows 
«/ as a function of 6/ for several values of Ucf for V/W = 
0.1 and n - 1.25. For a large Ucf, we recognize the three 
regimes mentioned above. A first-order phase transition oc- 
curs from the intermediate-valent regime to the «/ ^ 
regime for Ucf/W > 1.24. We observe hysteresis by increas- 
ing and decreasing e/ across the first-order phase transition 
point, and here we show the values of the state that has the 
lower energy. Figure 1(b) shows the valence susceptibility 
Xv = -drif/def as a function of ey. The valence susceptibility 
enhances around the boundaries of three regimes for a large 
Ucf - For a small Ucf, such a boundary is not clear and;^^v has 
a broad peak. Figure 1(c) shows the mass enhancement fac- 
tor l/Anikp) as a function of €f. In addition to the enhance- 
ment for Hf — > 1 as in the ordinary PAM, we find another 
region, that is, the intermediate-valent regime rif ^ 2 - n, in 
which the mass enhancement factor becomes large. This en- 
hancement, particularly, a peak as a function of e/, is not ex- 
pected for the PAM without Ucf. The large eff'ective mass in 
the intermediate-valent compounds Qr-YbAlB4 andj8-YbAlB4 
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Fig. 1. (Color online) e/ dependences of (a) n/, (b) x\, and (c) l/An{kp) 
for V/W = 0.1 and n = 1.25. Ucf/W = (dotted lines), 1 (dash-dotted lines), 
2 (dashed lines), and 3 (solid lines). 
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Fig. 2. (Color online) i/An(kj:) as a function of iif for V/W = 0.1 and 
n = 1.25. Ucf/W = (dotted line), 1 (dash-dotted line), 2 (dashed line), and 
3 (solid line). The thin Une is (1 
nf = 2 - n. 
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and the nonmonotonic variation in the effective mass under 
pressm^e in CeCu2Si2 may be explained by the present theory. 

To clearly observe the effect of Ucf on the mass enhance- 
ment, we show l/Anikp) as a function of «/ in Fig. 2. The 
thin line, which is almost overlapping with the Ucf - data, 
represents the mass enhancement factor (1 - n//2)/(l - «/) 
obtained for the PAM with Uc/ - and g = 1 . By increasing 
Ucf, 1 /An{kp) becomes large, particularly in the intermediate- 
valent regime ^ 2 - n. 



Fig. 3. (Color online) ;fv as a function of £f and Ucf with V/W = 0.1 for 
(a) n = 1.25, (b) n = 1.50, and (c) n = 1.75. The solid lines represent the 
first-order valence transition lines. The solid circles denote the critical points 
of the valence transition. The dotted lines indicate crossover lines determined 
by comparing the energies of the three extreme states (see text). 



In Fig. 3, we show the valence susceptibility ;t'v as a func- 
tion of 6/ and Ucf for n - 1.25, 1.50, and 1.75. In this fig- 
ure, we also draw the first-order valence transition lines and 
their critical points. The crossover lines, represented by the 
dotted lines, are determined by comparing the energies of the 
three extreme states: nf - 0, «y = 1, and «/ + «c/2 = 1 
with d = 0. The region where xv becomes large is captured 
well by the crossover hnes obtained by such a simple con- 
sideration. For n = 1.25, the first-order valence transition 
occurs only from the intermediate-valent regime to «/ ^ 
regime, while for n - 1 .75 it occurs only between the Kondo 
and intermediate-valent regimes, within the Ucf range pre- 
sented here. «/ in the intermediate-valent regime differs be- 
tween these two cases: «/ - 0.75 for n = 1.25 and «/ - 0.25 
for n - 1.75. The first-order transition seems to occur easily 
between very different states, that is, a crossover accompany- 
ing a large valence change tends to become a first-order phase 
transition. Between these two cases, for n - 1.50, both the 
transitions take place for Ucf/W > 2.88. Note that, since only 
the n = 1 .75 case is well investigated in previous studies, 
the first-order transition between the intermediate valent and 
«/ ^ regimes has not been elucidated. 
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Fig. 4. (Color online) 1 IMk^) as a function of e/ and Ucf with V/ W = 0. 1 
for (a) n = 1.25, (b) n = 1.50, and (c) n = 1.75. The lines and circles are the 
same as those in Fig. 3. 



Figure 4 shows the mass enhancement factor l/An(A:F) as 
a function of e/ and Ucf- A large mass enhancement occurs 
in the intermediate-valent regime in addition to the Kondo 
regime. Here, we note that the large mass enhancement oc- 
curs in the middle of the intermediate-valent regime. Thus, 
this enhancement is not due to the valence fluctuations. In 
CeCuiSia, the effective mass has a peak before the super- 
conducting transition temperature becomes maximum under 
pressure. If the system is in the Kondo regime at ambient 
pressure, passes the intermediate-valent regime under pres- 



sure, and finally reaches near a critical point, it is consistent 
with our theory provided the pairing interaction of supercon- 
ductivity is mediated by the valence fluctuations. Such a sit- 
uation can be realized, for example, for n - 1 .5 as is shown 
in Fig. 4(b). A similar discussion may also be applicable to 
CeCu2Ge2 if we can subtract the contributions of the mag- 
netic fluctuations. 

In summary, we have studied the extended periodic Ander- 
son model with Ucf by Gutzwiller approximation. We have 
found that three regimes, that is, the «/ ^ 0, intermediate 
valent, and Kondo regimes, are clearly defined for a large 
Ucf. Then, we have found that, in the intermediate-valent 
regime, the eff'ective mass is enhanced substantially. Accord- 
ing to the present theory, the large mass enhancement in the 
intermediate-valent regime indicates a large Ucf. Thus, our 

theory provides helpful information for searching a supercon- 
ductor with valence-fluctuation-mediated pairing. 



1) F. Steglich, J. Aarts, C. D. Bredl, W. Lieke, D. Meschede, W. Franz, and 
H. Schafer: Phys. Rev. Lett. 43 (1979) 1892. 

2) B. Bellarbi, A. Benoit, D. Jaccard, J. M. Mignot, and H. F. Braun: Phys. 
Rev. B 30 (1984) 1182. 

3) E. Vargoz and D. Jaccard: J. Magn. Magn. Mater. 177-181 (1998) 294 . 

4) H. Q. Yuan, F. M. Grosche, M. Deppe, C. Geibel, G. Spam, and 
F Steglich: Science 302 (2003) 2104. 

5) K. Miyake, O. Narikiyo, and Y. Onishi: Physica B 259-261 (1999) 676. 

6) Y. Onishi and K. Miyake: J. Phys. Soc. Jpn. 69 (2000) 3955. 

7) D. Jaccard, H. Wilhelm, K. Alami-Yadii, and E. Vargoz: Physica B 
259-261 (1999) 1. 

8) A. T. Holmes, D. Jaccard, and K. Miyake: Phys. Rev B 69 (2004) 
024508. 

9) T. M. Rice and K. Ueda: Phys. Rev. B 34 (1986) 6420. 

10) R Fazekas and B. H. Brandow: Phys. Scr. 36 (1987) 809. 

11) R. T. Macaluso, S. Nakatsuji, K. Kuga, E. L. Thomas, Y. Machida, 
Y. Maeno, Z. Fisk, and J. Y. Chan: Chem. Mater 19 (2007) 1918. 

12) S. Nakatsuji, K. Kuga, Y. Machida, T. Tayama, T. Sakakibara, 
Y. Karaki, H. Ishimoto, S. Yonezawa, Y. Maeno, E. Pearson, G. G. Lon- 
zarich, L. Balicas, H. Lee, and Z. Fisk: Nat. Phys. 4 (2008) 603. 

13) M. Okawa, M. Matsunami, K. Ishizaka, R. Eguchi, M. Taguchi, 
A. Chainani, Y. Takata, M. Yabashi, K. Tamasaku, Y. Nishino, 
T. Ishikawa, K. Kuga, N. Horie, S. Nakatsuji, and S. Shin: Phys. Rev. 
Lett. 104 (2010)247201. 

14) Y. Onishi and K. Miyake: Physica B 281-282 (2000) 191. 

15) S. Watanabe, M. Imada, and K. Miyake: J. Phys. Soc. Jpn. 75 (2006) 
043710. 

16) Y. Saiga, T. Sugibayashi, and D. S. Hirashima: J. Phys. Soc. Jpn. 77 
(2008) 114710. 

17) C. E. T. Gon^alves da Silva and L. M. Falicov: Sohd State Commun. 17 
(1975) 1521 . 

18) M. C. Gutzwiller: Phys. Rev 137 (1965) A1726. 



4 



